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We generalize the Nozieres-Schmitt-Rink (NSR) method to study the repulsive Fermi gas in 
the absence of molecule formation, i.e., in the so-called "upper branch". We find that the system 
remains stable except close to resonance at sufficiently low temperatures. With increasing scattering 
length, the energy density of the system attains a maximum at a positive scattering length before 
resonance. This is shown to arise from Pauli blocking which causes the bound states of fermion pairs 
of different momenta to disappear at different scattering lengths. At the point of maximum energy, 
the compressibility of the system is substantially reduced, leading to a sizable uniform density core 
in a trapped gas. The change in spin susceptibility with increasing scattering length is moderate 
and does not indicate any magnetic instability. These features should also manifest in Fermi gases 
with unequal masses and/or spin populations. 
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Since the early days of quantum many-body theory, 
the Fermi gas with a repulsive short range interaction 
has been used as the primary example of a Fermi liq- 
uid state[T]. The discovery of BEC-BCS crossover[2], 
however, shows that the ground state of this system 
is a molecular condensate, i. e., the Fermi liquid state 
is metastable. In the last two years, after the ref. [3] 
suggested the evidence of Stoner ferromagnetism, there 
has been increased interest in the nature of uncondensed 
Fermi gas (free of molecules) in the strongly interacting 
regime. Such systems have been referred to as the "up- 
per branch" Fermi gas, and the molecular condensate as 
the "lower branch" . 

Theoretical studies have found both ferromagnetic 
transition as well as the absence of it. [4] Though seldom 
emphasized, the upper branch Fermi gas in the strongly 
interacting regime has been studied by many experimen- 
tal groups[5, 6 at higher temperatures with different den- 
sities and trap depths. The key features in atom loss 
and energy maximum reported in ref. [3] also appeared 
in these earlier experiments. A puzzling feature is the 
presence of a range of scattering length (a s ) where the 
energy derivative is negative (d£/d(~aj 1 ) < 0) in ap- 
parent violation of the adiabatic relation of Tan[7]. Since 
the Fermi gas is unlikely to be ferromagnetic in the tem- 
perature regime of these earlier experiments, it leads to 
a natural and intriguing question on the nature of the 
repulsive gas in the strongly interacting regime. 

The key obstacle in theoretical studies of the upper 
branch Fermi gas is to find a proper mathematical de- 
scription of the "upper branch" . There is no precise for- 
mulation of it to the best of our knowledge. Fortunately, 
the meaning of upper branch is well defined in the high 
temperature regime, as the second virial coefficient 62 is 
made up of a bound state contribution and an extended 
(or scattering) state contribution, 6 2 = &| + ^2° ■ The 
upper branch corresponds to excluding the Hilbert space 



of molecules by setting Wf = 0. In addition, any descrip- 
tion of the upper branch Fermi gas must also recover the 
results 1 in the weakly interacting limit, as the energies 
of the molecules are very far below the continuum and 
hence can be ignored. 

Here we generalize the approach of NSR[S], which we 
call the excluded molecular pole approximation (EMPA) , 
to study the upper branch Fermi gas. It amounts to 
excluding the Hilbert space of molecules in a Guas- 
sian fluctuation theory [5], and obtaining thermodynam- 
ics within this truncated space. This approach recovers 
rigorously both the high temperature results and the re- 
sults of Galitskii in the weak coupling limit. Applying 
this method to lower temperature and strongly interact- 
ing regime, we find the following: (I) On approaching 
the resonance from the repulsive side at a fixed temper- 
ature T, the energy density £ attains a maximum at a 
positive scattering length (a sm ) prior to resonance, as 
seen in experiments. [3l |6] The theory also explains the 
subsequent fall in the energy density with increasing a s 
(violation of the adiabatic theorem of Tan). (II) The 
compressibility k attains a minimum at a sm (where £ is 
maximum). The small compressibility implies a core of 
almost uniform density at the centre of the trap. (Ill) 
The spin susceptibility x attains a maximum at the loca- 
tion of the energy maximum, i. e., at a sm ; it shows only 
a moderate variation over the entire range of a s , without 
any divergence indicative of a magnetic instability. 

EMPA for the Upper Branch Fermi Gas: Let us first 
recall that at low fugacity regime [ID] , the equation of 
state is n(T, /j,) = n (T, + dAP/dfi, where n (T, fi) is 
the density of an ideal gas, AP(T,/i) = T(^2/Xfz 2 h 2 is 

the interaction contribution to the pressure, A = yj 

is the thermal wavelength, and z — e M ^ T is the fugac- 
ity (h = ks = 1), fJ- is the chemical potential. The 
second virial coefficient 62 is made up of a bound state 
contribution b^ 1 and a scattering state contribution fo| c , 
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and — |.E&| = — (ma 2 ) -1 is the energy of the bound state, 
and rj is the phase shift. The interaction contribution 
to equation of state An(T, //) = n(T,fx) — n (T,[x) can 
therefore be written as An = An bd + An sc , where 



An Q (T, u) 




a = fed, sc (2) 



Next we recall that in the NSR approach[8], the inter- 
action contribution to the density An(/j,T) = n(T, fi) — 
n (T, fi) is 
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where ft is the volume, n (T,n) = 2^ fe n i? (£fc), is the 
density of a two-component ideal Fermi gas, n F (uj) = 
l/(e w / T + 1), a = £fe - e fc = /c 2 /2m, n B (w) = 
l/(e w / T — 1), and M(uj + ,q) is negative inverse of the 
two particle T-matrix in the medium, of the form 

J L(uj + ,q) (4) 
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j(k;q) = l — n F (t;z +k ) — n F (l;z_ k ) describes Pauli block- 
ing of pair fluctuations. In the extreme dilute limit, 
j(k;q) reduces to 1, — M^ 1 to the inverse two-body T- 
matrix, and the phase angle £(to,q) = a,rgM(u> + ,q) to 
the negative of the two-body phase shift rj(ui — w(q)), 
where uj(q) = q 2 /4m — 2//. 

For a given q, the value of ((ui,q) depends on the 
location of branch cut and poles of M~ 1 (u) + ,q). It is 
clear from Eq.Q and ([5| that the branch cut is given 
by uj > Lu(q). Should M~ 1 (uj + ,q) have a pole, say at 
u>b(q) < w (<?)> then we have 



ImL(w, q) 



u > C(^,<?) = tan 1 

uj b (q) <uj < u(q), ((uj, q) = -7T, 

w<w 6 (g), c( w ,g) = o. 

Otherwise, Eq.([7]) and Q are replaced by 
u><u(q), ((u,q) = 0. 
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Eq.([3} can then be written as An(T, //) = An bd (T, 
An sc (T,u), 

An M (r )A i) = -ix;»BMg))^ (10) 
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FIG. 1. (Color online) (a) Position of pole ujb(q) in cj — q 
plane for a given a s : The solid curve denotes the curve 
cu(q) = q 2 /(4m) — These pole position uJb(q) is the so- 
lution of the equation ReM(cj + ,q) = (see eqn. Q), where 
Pauli blocking is described by y(k; q). For a s > 0, the matrix 
M (lo + ,q) of a two-body system will have a pole of energy 
— \Eb\ below Lj(q), (dotted line). In a many body system, 
Pauli blocking will suppress formation of molecular bound 
states. The suppression is strongest for pairs with total mo- 
mentum q = and is less strong for larger q. As a result the 
pole position changes to those indicated by the dashed blue 
curve, (b) The red curve is the critical scattering length a c s (q) 
at T = 3E F . (See (B) in Summary of Results). For a given 
a s , a fermion pair with total momentum q (referred simply 
as "g-pair") can have a bound state only when a s < <is(q), 
i.e., to the left of the red curve. As a s increases, such that 
a s crosses a c 3 (q) from left to right, a q-pair will lose its bound 
state, and the energy of the scattering state of this pair will 
jump downward abruptly (see fig. EJ. 



That we use the same superscript in Eq.( 10 ) and ( 11 ) as 



in the high temperature case is because they reduce to 
Eq.|2]) in the low fugacity regime. Thus, by continuity, 
the extension of the upper branch Fermi gas to lower tem- 
perature is to exclude the contribution from the molecu- 
lar bound pole term ( Eq.(10)) from An(T, //). Hence the 



name EMPA. The equation of state within EMPA is then 



n(T, M )=n (T, M )+An sc (T,/i). 



(12) 



Considering a system with fixed density n (which defines 
Ep and kp) and inverting the relation n = n(T,/i) to 
obtain /1 = /x(n,T), one can obtain all thermodynamic 
potentials as a function of n and T. [IT] It is also straight- 
forward to describe spin polarized states in the upper 
branch by including a magnetic field h. 

Note that Eq.(Jl2| involves only integrating over the 
area uj > uj(q) (i. e., above the solid curve in Fig. [l|a)) 
with an integrand given explicitly by Eq.([6|. There is 
no need to obtain the pole structure as far as evaluating 
Eq.(12) is concerned. There is, however, a close connec- 



tion between the interaction energy of scattering state 
and the presence of a pole. Understanding the distribu- 
tion of poles in the ui-q plane is, therefore, essential for 
the elucidation of the results to be presented below. 
Summary of Results: 

(A) Phase diagram: Fig. [2] displays the "phase diagram" 
of the upper-branch Fermi gas. All regions except the 
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FIG. 2. (Color online) Upper branch "phase diagram". The 
point K corresponds to (— l/k F a s = — 0.435, T = 2.85E F ). 
The solid blue line ending at K is a locus of states with a 
vanishing compressibility. The dashed curve starting at K 
shows a sm (see text). The region hatched in dark grey is 
mechanically unstable. Tan's adiabatic theorem is violated in 
the region shaded in light grey. 

region in dark grey are stable (k,X > 0)- The solid line 
that ends at K is the boundary of vanishing compress- 
ibility k — 0. The dashed line above K describes a sm 
where the energy attains a maximum at a fixed temper- 
ature. Across this line, fi, k, x, and energy density £ 
are continuous but their slope undergoes sharp changes. 
These discontinuous slopes, however, may disappear if 
beyond Gaussian fluctuations are included. Crossing the 
solid line below the point K, the quantitites /i, P, k, 
X, and £ undergo discontinuous changes; the system 
is mechanically unstable. The white and light grey re- 
gions correspond to regimes with d£ /d(— aj 1 ) > and 
d£/d(— aj 1 ) < respectively. In the light grey region, 
Tan's adiabatic theorem is not applicable (see below). 
(B) Energy density £ : Fig. [3]ja) shows the behavior of 
energy density £ as a function of kp a s at T — 3Tp. It 
exhibits a maximum at kpa sm = 2.61, (which falls on the 
dashed line in Fig. [2]). Such a maximum feature is consis- 
tent with the early observation by Salomon's groupjB] at 
high temperatures, as well as in ref.[3] at lower temper- 
atures. The maximum behavior implies that there is a 
region of kpa s (the light grey region in Fig. [2]) where the 
the adiabatic theorem, d£ /d{— aj 1 ) > is violated. The 
resolution is that the relation between d£/d(— aj 1 ) of 
the scattering state and the contact density is ill-defined 
at the scattering length where a molecular bound state 
disappears. 

This is best seen in the two-body case (see Fig. [2|, 
where the energy of the scattering state of a fermion pair 
with total momentum q (referred to as "q-pairs") jumps 
suddenly downward when of a s passes a critical value 
(( a s) _1 = in this case) at which the molecular bound 
state on the side a s < (a^) -1 disappears. In the many- 
body case, due to Pauli blocking, (j(k; q) ^ 1), different 
q-pairs will form bound states at different critical scat- 
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FIG. 3. (Color online) (a) Energy density (£) and (b) com- 
pressibility («;) and susceptibility (x) as a function of the scat- 
tering length at T = 3E F . All quantities are measured in 
units of their respective values of the non-interacting gas (in- 
dicated by the subscript o) at the same temperature. 

tering length a° s {q) (which is the lowest value of a s such 
that the equation ReM(w, q) — has a solution). Since 
Pauli blocking effect is strongest for the q = molec- 
ular bound state (Fig. [ija)), and is less significant as 
q increases, (a^q)) -1 is largest at q = and decreases 
monotonically as q increases. The behavior of a c s is shown 
in Fig. [TJb), and a sm = a c s {q = 0). 

That d£ /d{— aj 1 ) < for upper branch Fermi gas suf- 
ficiently close to resonance is now clear. As a s passes 
through a c s (q) from below, the molecular bound state of 
a q-pair disappears because of Pauli blocking. Up on this 
disappearance, the energy of the scattering states of this 
pair suddenly jumps down, thereby causing the energy 
to decrease. As a s continues to increase, q-pairs with 
successively higher total momentum q lose their bound 
states, inducing a successive downward jump in the en- 
ergies of the scattering states of these pairs, and hence 
a negative derivative d£/d(— a^ 1 ) < 0. Note that since 
a sm is determined only by Pauli blocking, it should be a 
universal function of T and n i.e., kpa sm = f(T/Ep), 
where / is a dimensionless function (dashed line in 
Fig.§. 

Our explanation above might lead one to think that 
the energy decreasing process will cease when no more 
q-pairs lose their bound states, which occurs at a s = oo. 
What actually happens, however, is that the minimum of 
£ as a 8 increases beyond a sm (which signifies the ceasing 
of energy decrease) occurs at a scattering length prior to 
resonance. The reason is that in order to have an energy 
decrease caused by the scattering state of a q-pair, this 
pair state has to be occupied. At lower temperatures, 
the probability of occupation of such pair states is low 
especially for those pairs with high q, thereby causing the 
energy decrease to cease at an (a s ) m i n prior to resonance. 
As T increases, ({a s ) m i n )~ 1 approaches 0. 
(C) Compressibility n: As a s increases, a repulsive Fermi 
gas is expected to become less compressible. For tern- 
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FIG. 4. (Color online) The discontinuous change of the energy 
of the scattering state (solid line) of a two body system up on 
the disappearance of the molecular bound state (double line). 
|10| Similar phenomena occur in a q-dependent fashion in the 
many body setting. 
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FIG. 5. (Color online) Comparison of densities of strongly 
interacting and weakly interacting gases in a spherical trap. 
Ep corresponds to the density at the trap centre. For each 
temperature, the number of atoms in both the strong and 
weak cases is the same. A "flat-top" density profile is evident 
in the strong case, and becomes more pronounced at lower 

/ 2E° 

temperatures. The radius r is in units of y (^f-) [ft>t - trap 
frequency) . 

peratures above that of point K in Fig. [2j k attains a 
minimum at a s = a sm (see Fig. [3](b)) . Our calculation 
shows, for temperatures lower than that of K, k — > as 
one approaches the solid line (Fig. [2]) from the left. The 
system behaves like a hard core Fermi gas with a core 
size close to inter-particle spacing. There is, however, 
an important difference between a hard core Fermi gas 
with core size equal to a s ss kp 1 and the actual atomic 
Fermi gas. In the former case, the effective range is also 
of order kp , whereas the effective range in atomic gases 
is much less than the inter-particle spacing, independent 
of the value of a s . The diminished compressibility has 
a dramatic effect on the density profile. This leads to 
clouds with little variation of density at the centre, an 
effect that becomes more pronounced at lower tempera- 
tures (see Fig. [5]) . 

(D) Spin susceptibility \ : Fig- [3j^b) also shows the spin 
susceptibility x at T = 3Tp. Note that \ changes 



by at most 40 percent for the entire kpa s range, and 
only moderately in the experimentally relevant range 
0.5 < kpa s < 2. We do not see a diverging suscepti- 
bility indicative of a magnetic transition. 

Although our discussions focused on the equal-mass 
spin-i Fermi gas, the properties enumerated arise mainly 
from Pauli blocking. These features should, therefore, be 
generic to other upper branch Fermi gases, such as those 
with unequal masses or unequal spin populations. 
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